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PLANE WAVE PROPAGATION IN A LINEAR, HOMOGENEOUS, ANISOTROPIC 
DIELECTRIC MEDIA ("CRYSTAL OPTICS"): 

Our objective here is to formulate a general approach to the subject of wave propagation in 
anisotropic dielectrics which makes use of ideas famihar from other branches of 
mathematical physics — viz. the "eigenvalue problem.". For reasons that will soon 
become abundantly clear, treatments of "crystal optics" focus on the behavior of the 
dielectric displacement vector, D(r, CO ) rather than on the electric field vector. For non- 
magnetic dielectrics the components of the dielectric displacement are usefully represented 
as the Cartesian coordinates coordinates of figure called the "ellipsoid of wave normals," 
the "optical indicatix," the "index ellipsoid" or the "reciprocal ellipsoid." 



Since the stored electrical energy is given by 

f/, = (l/2)DE=(l/2)D-8"'-D 

we can write in the principle axis system 
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where {e^ , £,,,£,} are principle axis values of the dielectric 
constant tensor. 




Kaiser S. Kunz in 1977 presented a similar treatment of this problem in a paper enti tiled "Treatment of optical 
propagation in crystals using projection dyadics," Am. J. Phys., Vol. 45, 1977, pp. 267-269. 

Perhaps the most authoritative treatment of "Crystal Optics" is found in Max Born and Emil Wolf, Principle of 

Optics, Pergamon Press (1986), Chapter 14. r. victor jones, December 19, 2000 
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Now let us first combine Equations [ I-8a ] and [ I-8b ] to obtain a generalized Helmholz 
equation for a homogeneous anisotropic dielectric 

Vx(vx[r'((o)-D(r,(o)]) = (oV« iXr,co) [V-1 ] 

In general, we would hope to be able to find a set of eigenmodes of the homogeneous 
problem that would satisfy the scalar eigenequation 

V'D"'\r,03) + nl^^kl D'"\r,a))=0 [V-2] 

where n^^^ is the effective index of refraction of the a-th eigenmode. We can be quite 

specific for the case of plane wave eigenmodes where we suppose that all fields have a 
spatial dependence exp(+ / k r) . From Equation [ I-8c ] we see that D(r, co) must, in 
general, be orthogonal to k so that we may write 

D (k, CO) = D^''(k, CO) t'-'Xk, co) + D'''(k, co) t^'^(k, co) [ V-3 ] 

where the t*°' (k)'s are polarization unit vectors which are orthogonal to k , the unit vector 
parallel to k . Since the components of 8 in the general case may be complex ~ e.g., in the 
case of magneto-optical media ~ the eigenmodes may be polarized along "complex 
directions" ~ e.g., "screw axes" defining the sense of circular polarization ~ and we must 
use great care in all vector manipulations. Thus, we define a set of adjoint or conjugate unit 
vectors by means of the relationships 

k-k=l; k-t('"(k)=0; t(<"(k)-k = 0; ^''^(k)- t<''''(k)= 5„„, [ V-4 ] 

Thus, with k = k k = k* k (so that k k = |k|^ ) the representation for D(k) in Equation 
[ V-3 ] automatically satisfies Equation. [ I-8c ]. A key problem is the representation of the 
ubiquitous vector operation 



R. Victor Jones, December 19, 2000 



On Classical Electromagnetic Fields Page 44 

3 (r)-F(r)= V x(v xF(r)) = W •F(r)-V'F(r) [ V-5 ] 

which appears in Equation [ V-1 ]. For plane wave representation of the operator becomes 

3(k) = (k-k)i-kk = |k|'{i-kk}. [V-6] 

Since 

1 = k k + t''\k) V\k) + t*" (k) t*" (k) [ V-7 ] 

the plane wave representation of the operator simplifies to 

3 (k) = |k|' {t^'\k) t^'\k) + t'"' (k) t'"' (k)} [ V-8] 

Using this representation of the 3 (k) operator and the representation for D(k) in Equation. 
[ V-1 ], the generalized Heknholz equation ~ i.e. Equation [ V-1 ] ~ can be written 

ikf {D<"(k) v\k) ■ r • t^'^(k)] + (k) (k) [t (k) ■ ■ t^'^(k)] + 

D (k) t^''( k) [ t *'\k) • 8"' • t (k)] + D^'' (k) t^'' (k) [t*'' (k) • e"' • t''' (k)]} [ V-9 ] 
= 0) Vo{D''\k) P'(k) + D*^> (k) F(k) } 

Crucial point ; To obtain an eigenvalue equation we need to choose the eigenvectors 
t^°^(k) so that 

t^'\k) • e~' • t^" (k) = t^'' (k) • r' • t*"(k) = [ V-10 ] 

If we can find eigenvectors defined in this way, then Equation [ 1-9 ] becomes an 
eigenvalue equation with eigenvalues - i.e., the inverse refractive indices - given by 

[n^(k)]"' = k^/|k|'=8o P\k)-r'-F\k) [ V-lla] 
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[n,(k)]" =kl/\kf = e,V'\k)-r''e\k) [V-llb] 

These results ~ i.e., Equations [ I-IO ] and [ I-ll ] ~ are a complete formal solution of the 
problem. However, they are difficult to apply in the general case and an addition 
relationship ~ viz., the Fresnel equation of wave normals ~ is found to be extremely useful 
as the starting point for actual computations. For plane waves. Equation [ V-1 ] can be 
rewritten as 

[|k|' i - (koVeo) e] -£(1^) = k k •E(k) [ V-12 ] 

Multiplying this equation through by k • pk|^ i - (kg/eo) 8 j • we obtain 

k £(1^)= {k- [ikf i -(kl/E,)eY ■ k} k- E(k) [ V-13a ] 

or l=k-[|kf i-(ko/eo)8]"'-k [V-13b] 

From this result we may develop two important relationships. Using the principal axes 
coordinates of the dielectric tensor, we can write 



[ikp 1 - (koVeo) e] = - (ko/eo) \ a, [ V-14a ] 

X= a 

and [|k|' 1 - (klle,) e]"' = f^(\kf - (k^ej E(,))"' a, [ V-14b ] 



Therefore, Equation [ I-14b ] becomes 



(k.aj(k.aj (k.a,)(k.aj (k.a)(k.aj 1 
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where nl = £(x)/^o ■ Since a^^ = 1 (or ^(k • )( '^)-^ ) we may also write 

X=a X=a 

Equation [ V-15a ] as 

(k-aj{k-aj {k-a,){k-a,) {k-aj{k.aj 



1 1 



1 1 



+ ■ 



1 1 



= 



[ V-15b] 
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This latter expression is the famous Fresnel equation of wave normals. 

APPLICATIONS OF THE FORMAL SOLUTION: 
Uniaxial Dielectric Crystals : 

For an optical material with uniaxial symmetry, the inverse dielectric tensor in the 

principal axes system must have the form 



[V-16] 



Thus, Equation [ I-19b ] becomes 
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cos 13 r = o 



[V-17] 



so that 



=nl\ and =n'l cos^i^ + nn"^ sin^-d 



[V-18] 



[^l H n 

2" In its commonly used form, the Fresnel equation becomes ~ + ~ ~[ ~ ~ ^ 

2 2 2 2 2 2 

n n n 

21 In this instance there is no need to trouble ourselves about conjugate unit vectors. 

R. Victor Jones, December 19, 2000 



On Classical Electromagnetic Fields 



Page 47 



where k- 3^ = 00813^. The subscript "o" identifies the "ordinary" mode and the 
subscript "e" the "extraordinary" mode. These results are usually plotted as follows: 




where the intersections of the k vector yield the "ordinary" and "extraordinary" 
velocities of propagation for a given k . 

Further, if we take 

k = sin(p sini3 + cos(p sini3 a ^ + cosi3 a^ [ V-19 
it is a bagatelle to show that 

t*°'(k)= r^'Ck) = cos(p a^ - sincp a, [ V-20a ] 

t*'' (k) = t (k) = cos-d (coscp a^ + sincp a J -siniJ a, [ V-20b ] 

and that these quations are consistent with Equations [ V-10 ] and [ V-11 ]. 
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See, for example. Section 82 in L. D. Landau and E. M. Lifshitz, Electrodynamics of Continuous Media, 
Pergamon Press (1960). 
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i+tY(a„a,-a,aj]8-^ [V-21] 



If we introduce the conjugate principal axes 

a.=a_=(l/^)(a„+ta,) 

5_ = a, =(l/^)(a,-taj [V-22] 

we obtain the dielectric dyadic in the so called normal form -- viz. 

= [(1 - Y ) a, a, + (1 +Y ) S_ a_ + 5„ S„ ] e"^ [ V-23 ] 

Again from Equation [ V-15b ] it is trivial to show that 

=[1±Y cosT^]8"' [V-24] 

Using the resolution of k as given in Equation [ V-4 ] we may show that 

P (k) = (^^ ^ exp(- 1 (p) [ 1 + cos 1^] a^ - exp(i (p) [1 ± cos a_ ± ^2 sini^ a,, } [ V-25a ] 

and 

t-(k)= [^^j{-exp(i(p)[l+cosi^] a+ +exp(-i(p)[l±cosi^] a_ + >72 sini^ a,,} [ V-25b ] 

ENERGY FLOW IN ANISOTROPIC MEDIA 

As previously noted, the content of Equations [ I-IO ] and [ I-ll ] represents in some 
sense a complete formal solution of the wave propagation problem. However, from a 
practical point of view it is essential to consider how energy propagates in anisotropic 
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media. To that end, we note that the time averaged Poynting vector associated with a 
given plane wave- viz. 



S (k, (o) = ^ E (k, (o) X H* (k, (o) 



[V-26] 



propagates in a direction s (k, (oj — conventionally designated the "ray" direction 
which is orthogonal to both E (k, co) and 11 (k, O)) as shown below. 




The time averaged total stored energy is given by 
1 r - /r 



t/(k, co) = - [e (k, (o) ■ D* (k, co) + B (k, co) ■ H* (k, co)] 



= ^k E(k,co)xH* (k,co) = ^'k S(k,(o) 



[ V-27 ] 
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and, thus, we see that the "ray" or "energy flow" velocity, v^^y , for a given s (k, co) i 
given by 



IS 



1 I ^ \ 1 ^ \ 
— = — k s k, CO = k s k, CO 



[ V-28 ] 

' ray ™ ' phase 

We write the time averaged Poynting vector associated with a given eigenmode as 



S"" (k) = i E"" (k) X H'-i* (k) = E"" (k) X [k X E"'* (k)' 



[ V-29 ] 



where E*"' (k) is the electric field associated with the eigenmode. Using Equations 
[ V-3 ], [ V-10 ] and [ V-U ] this field can be expressed as 

E<"' (k) = {[t'"\k) • 8"' • t<"\k)] V'ik) + [k-r' ■ t'"\k)] k} D'"'(k) 



r ^ ^.2 {t^'^kk) + Eo [«(a)(k)]'[k- e"' • t'^^k)] k} D'^^k) 
[«(c)(k)J 



[ V-30 ] 



^ {t'">(k) + g^">(k)k}D'°\k) 



where g<"\k)= eJ.^.^Ck)] [k.rM<^>(k)] = ^1^^^^^^ 
parameterization, the modal Pojmting vector can be expressed as 



Using this 



[V-31] 
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and the associated ray vector as 

k-t<'"(k)^^'"(k) k-P'(k)^^'''(k) 



Jk-t<''>(k)g^<'>(k)| Ji+ig^^'ck)!' 



If we take k and t'^^k) as reference directions, the ray vector, s*"' (k), lies in the 
plane containing k and t'^'''(k) at an angle T|*'"(k)= -tan'^fg^^^k)] with respect to the 
direction k. 

For an optical material with uniaxial symmetry, we may use Equations [ V-16 ], 
[ V-19 ] and [ V-20 ] to evaluate g'''^(k). In particular, we may easily see that for the 
ordinary mode 

tami*°\k) = -g^°'(k) = [ V-33a ] 

or s^°' (k) = k and for the extraordinary mode 



sin-dcosi3^ (£± 

e;^' cos^-d + e,;^ sin^iJ' 



tanri^^^k) =-/\k) = ^y.]:''^,^^^ . [ V-33b ] 
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